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Abstract: We present measurements of radioactive contamination in the high-resistivity silicon 
charge-coupled devices (CCDs) used by the DAMIC experiment to search for dark matter particles. 
Novel analysis methods, which exploit the unique spatial resolution of CCDs, were developed to 
identify a and j3 particles. Uranium and thorium contamination in the CCD bulk was measured 
through a spectroscopy, with an upper limit on the (^^^Th) decay rate of 5 (15) kg^^ d^' at 
95% CL. We also searched for pairs of spatially correlated electron tracks separated in time by 
up to tens of days, as expected from ^^Si or -^*°Bi sequences of jS decays. The decay 
rate of ^^Si was found to be kg^^ d^' (95% Cl). An upper limit of ~35 kg^* d^^ (95% CL) 
on the ^'^^Pb decay rate was obtained independently by a spectroscopy and the j8 decay sequence 
search. These levels of radioactive contamination are sufficiently low for the successful operation 
of CCDs in the forthcoming lOOg DAMIC detector. 

Keywords: Dark Matter detectors (WIMPs, axions, etc.); Solid state detectors; Very low-energy 
charged particle detectors; Search for radioactive and fissile materials. 
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1. Introduction 

The DAMIC (Dark Matter in CCDs) experiment []T], 0 employs the bulk silicon of scientific-grade 
charge-coupled devices (CCDs) fo defecf coherenf elasfic scattering of Weakly-Inferacfing Massive 
Particles (WIMPs) — pufafive yef-fo-be-discovered particles which may explain fhe dark matter in 
fhe universe 0-i- By virtue of fhe low readouf noise of fhe CCDs and fhe relatively low mass 
of fhe silicon nucleus, DAMIC is particularly sensitive fo low mass (<20GeV/c^) WIMPs, which 
induce nuclear recoils of keV-scale energies. 

As for any direcf dark matter search, fhe ulfimafe sensifivify of fhe experimenf is defermined by 
fhe rale of radioactive background fhal mimics fhe nuclear recoil signal from WIMPs. Al fhe lowesf 
energies (<lkeV), CCDs, as well as any olher presenl defector fechnology, lack fhe capabilify 
of effective discriminafion befween signals produced by nuclear recoils and fhose from ionizing 
eleclrons, making energy deposifs in fhe active fargef from j8s and y-rays pofenlial backgrounds for 
a WIMP search. 

To suppress such pofenlial backgrounds a breadfh of slrafegies have been adopfed. Direcf dark 
matter search experimenls are deployed in deep underground laborafories, fo eliminafe cosmogenic 
backgrounds. DAMIC is located in fhe SNOLAB laboratory, 2 km below fhe surface in fhe Vale 
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Creighton Mine near Sudbury, Ontario, Canada. Dedicated screening and selection of detector 
shielding materials, as well as radon-suppression methods, are extensively employed to decrease 
the background from radioactive decays in the surrounding environment. However, it is the ra¬ 
dioactive contamination of the active target that will often determine the feasibility and scalability 
of different technologies for rare-event searches. Thus, the measurement of the intrinsic contami¬ 
nation of the detector is fundamental, especially for solid-state devices, whose active target cannot 
be purified after fabrication. Particularly relevant for DAMIC and other silicon-based experiments 
(e.g. SuperCDMS [^) is the cosmogenic isotope ^^Si, which could be present in the bulk of the 
detector. Its j8 decay spectrum extends to the lowest energies and may ultimately become an iiTe- 
ducible background. 

The dark matter search will be performed with DAMIClOO, a detector with 100 g of sensitive 
mass. We aim for background rates in the WIMP search region of ~2 keV^^ kg^'d^\ comparable 
to those of state-of-the-art ionization detectors [0, This will allow us to probe, after one year 
of operation, WIMP-nucleon spin-independent interaction cross-sections as small as 10^^ pb for 
WIMPs with masses as low as 2 GeV/c^. To achieve this goal, the total decay rate of j3 emitters in 
the bulk of the CCD must be ~300 kg^'d^ ’ or lower. 

In this paper, we present novel techniques for the measurement of radioactive contamination 
in the bulk silicon and on the surface of DAMIC CCDs. We exploit the superb spatial resolution 
of the CCD to derive distinctive signatures for a and j3 particles, and to identify radioactive decay 
sequences with time separation of up to weeks. 

We describe the properties of the CCDs and the DAMIC setup at SNOLAB in Section ^ Char¬ 
acteristic features of a-induced signals, and limits on uranium and thorium contamination through 
a spectroscopy are presented in Section 0. We detail in Section || a search for spatially correlated 
P decay sequences and corresponding results on ^^Si and ^'®Pb contamination. Conclusions are 
drawn in Section 0. 

2. The DAMIC detector 

2.1 Characteristics of DAMIC CCDs 

The DAMIC CCDs were developed at Lawrence Berkeley National Laboratory Microsystems 
Lab 10, starting from an existing design for the Dark Energy Survey (DES) camera (DECam) m- 
They feature a three-phase polysilicon gate structure with a buried p-channel. The pixel size is 
15 /tm X 15 /tm and the active region of the detector is high-resistivity (10-20 kflcm) n-type sil¬ 
icon hundreds of micrometers thick. The high-resistivity of the silicon allows for a low donor 
density in the substrate (~10*^ cm^^), which leads to fully depleted operation at reasonably low 
values of the applied bias voltage (~40V for a 675 /rm-thick CCD). The CCDs are typically 8 
or 16 Mpixels, with surface areas of tens of square centimeters. Eigure |I] shows a cross-sectional 
diagram of a CCD pixel, together with an illustration of the WIMP detection principle. 

When operated at full depletion, ionization produced in the active region will be drifted along 
the direction of the electric field (z-axis). The holes (charge carriers) will be collected and held near 
the p-n junction, less than 1 /im below the gates. Due to the mobility of the charge carriers, the 
ionized charge will diffuse as it is drifted, with a spatial variance that is proportional to the carrier 


- 2 - 


3-phase 
CCD structure 

Poly gate 
electrodes 


(10 kn-cm) 




Photo¬ 
sensitive 
volume 
(500 pm) 
1 


—Buried 
p channel 




^ / 
Transparent 
rear window 


_ Bias 
voltage 


coherent elastic scattering 



(a) A CCD pixel 


(b) WIMP detection in a CCD 


Figure 1. a) Cross-sectional diagram of a 15 /dm x 15 jim pixel in a fully depleted, back-illuminated CCD. 
The thickness of the gate structure and the backside ohmic contact are <2 jum. The transparent rear window, 
essential for astronomy applications, has been eliminated in the DAMIC CCDs, b) Dark matter detection in 
a CCD. A WIMP scatters with a silicon nucleus producing ionization in the CCD bulk. The charge carriers 
are then drifted along the z-direction and collected at the CCD gates. 



(a) Portion of a DAMIC image 


(b) Emission of Si fluorescence X-ray 


Figure 2. a) 50x50 pixel portion of a CCD image, taken when the detector was at ground level. Different 
kinds of particles are recognizable (see text). For better contrast, only pixels with deposited energy >0. 1 keV 
are represented in color, b) Event with two nearby clusters detected after illuminating the CCD with a ^^Fe 
source. The 1.7 keV cluster is a photoelectron (pe) from the absorption of a Si fluorescence X-ray, emitted 
following photoelectric absorption of the incident 5.9 keV Mn X-ray in a nearby site. 


transit time. Charge produced by interactions closer to the back of the CCD will have longer transit 
times, leading to greater lateral diffusion. The lateral spread (width) of the charge recorded on 
the CCD x-y plane may be used to reconstruct the z-coordinate of a point-like interaction [@]. For 
extended tracks, e.g. from electrons and muons, this effect leads to a greater width when the track 
is closer to the backside, which provides information on the track orientation. 
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The ionized charge is held at the CCD gates until the charge is read out. During readout, 
the charge is transferred vertically from pixel to pixel along each column by appropriate clock¬ 
ing of the 3-phase gates (“parallel clocks”), while higher frequency clocks (“serial clocks”) move 
the charge of the last row horizontally to a charge-to-voltage amplifier (“output node”). The in¬ 
efficiency of charge fransfer from pixel fo pixel is as low as 10^^ and fhe readouf noise for fhe 
charge collecfed in a pixel is |Q]. Since on average 3.6 eV is required fo ionize an elecfron 

in silicon, fhe readouf noise corresponds fo an uncerfainfy of ~7 eV in deposifed energy. The im¬ 
age can then be reconstructed from the order in which the pixels were read out, and contains a 
two-dimensional stacked history (projected on the x-y plane) of all particle interactions thr'ough- 
out the exposure. For rare-event searches, it is advantageous to take long exposures (~8 hours in 
DAMIC) in order to minimize the number of readouts, and thus the number of pixels above a given 
threshold due to noise fluctuations. Note that even with these long exposures the CCD dark current 
(<0.1 e^pix^May^^ at the operating temperature of ~140 K) contributes negligibly to the noise. 

Figure [2(a)| shows examples of particle tracks as recorded by a DAMIC CCD. Low energy 
electrons and nuclear recoils, whose physical track length is <15 /rm, produce “diffusion limited” 
clusters, where the spatial extension of the cluster is dominated by charge diffusion. Higher energy 
electrons, from either Compton scattering or j3 decay, lead to extended tracks. Alpha particles in 
the bulk or from the back of the CCD produce large round structures due to the plasma effect El 
(Section ^71. Cosmic muons pierce through the CCD, leaving a straight track of minimum ionizing 
energy. The orientation of a muon track is immediately evident from its width, since the end¬ 
point of the track on the back of the CCD is much wider than the end-point at the front due to 
charge diffusion. In Figure |2(b)L we show a Mn X-ray interaction from a ^^Fe source which further 
demonstrates the superb spatial resolution of the CCD. Usually, the full 5.9 keV energy of the Mn 
Ka X-ray is deposited as a single cluster. Rarely, a fluorescence X-ray is emitted following phofo- 
elecfric absorption of the primary X-ray. The fluorescence X-ray may fravel far enough wifhin fhe 
CCD (a few attenuation lengths) to deposit its energy away from the first X-ray interaction, leading 
to two separate clusters. 

The detectors present an excellent linearity and energy resolution (55 eV RMS at 5.9 keV) |@] 
for electron-induced ionization, as measured with X-ray sources. The CCD response to a parti¬ 
cles was calibrated with an ^‘^'Am source; the energy scales from a-induced and electron-induced 
ionization were found to be within 3%, with an a energy resolution of 50keV at 5.5 MeV EUl- 

2.2 Setup at SNOLAB 

Most of the infrastructure for the DAMIC dark matter search is already installed in SNOLAB 
(Figure A packaged CCD (2kx4k, 8 Mpixel, 500 /tm-thick) is shown in Figure 0(a). The 
device is epoxied to a high-purity silicon support piece. The CCD clocks and output node signal 
travel on a Kapton cable, appropriately shaped for wire bonding. The cable is also glued to the 
silicon support. A copper bar facilitates the handling of the packaged CCD and its insertion into a 
slot of an electropolished copper box (Figure 0(b)). The box is cooled to ~140 K inside a copper 
vacuum vessel (~10^®mbar). An 18cm-thick lead block hanging from the vessel-flange shields 
fhe CCDs from radiation produced by the electronics card (Vacuum Interface Board, VIB), also 
located inside the vessel (Figure 0(c)). The CCDs are connected to the VIB through the Kapton 
flex cables, which run along the side of the lead block. The processed signals then proceed to the 
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Figure 3. a) A packaged DAMIC CCD. b) The copper box housing the CCDs, c) Components of the DAMIC 
setup, ready to be inserted in the vacuum vessel, d) The vessel inside the lead castle, during installation of 
the polyethylene shield. 


data acquisition electronic boards. The vacuum vessel is inserted in a lead castle (Figure 0(d)) 
which shields the CCDs from ambient 7 -rays through at least 21 cm of lead. The innermost inch of 
lead comes from an ancient Spanish galleon and has negligible ^^®Pb content, strongly suppressing 
the background from bremsstrahlung 7 s produced by ^^°Bi decays in the outer lead shield. A 
42 cm-thick polyethylene shielding is used to moderate and absorb environmental neutrons. 

Three 2kx4k, 8 Mpixel CCDs were deployed in February 2014 as part of DAMIC R&D ef¬ 
forts. Two of them were originally designed for astronomy, with a transparent indium-tin-oxide 
(ITO) coating deposited on the backside after thinning the CCD to 500 /^m . The third CCD was 
optimized for DAMIC by maximizing its mass (the CCD is un-thinned, 675 /im-thick) and mini¬ 
mizing radioactive contamination (the ITO layer containing j3-radioactive ’^^In is eliminated). The 
500-;Um CCDs are inserted in adjacent slots of the copper box, with copper plates above and below. 
The 675-/xm CCD is in a lower slot of the box, separated from the other CCDs by ~1 cm of copper. 
DAMIClOO will consist of 18 CCDs, each of 4kx4k pixels and 675 /xm thickness. 

2.3 CCD image reduction and data samples 

Clusters of energy deposits are found in the acquired images with the following procedure. First, 
the pedestal of each pixel is calculated as its median value over the set of images. The pedestals are 
then subtracted from every pixel value in all images. Hot pixels or defects are identified as recurrent 
patterns over many images, and eliminated (“masked”) from the analysis (>95% of the pixels were 
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Table 1. Analyzed data sets. The background rate refers to electron-like tracks in the 675 /rm-thick CCD. 


Start date 

End date 

Live-time 

day 

Pixel integration time 
/IS 

Background rate 

2014/06/06 

2014/07/07 

28.7 

0.6 

45.9±0.7 

2014/07/11 

2014/08/21 

36.1 

40 

5.1±0.2 

2014/09/03 

2014/09/29 

20.7 

40 

4.8±0.3 


deemed good). Pixel clusters are selected as any group of adjacent pixels with signals greater 
than four times the RMS of the white noise in the image. The resulting clusters are considered 
candidates for particle interactions. Relevant variables (e.g. the total energy by summing over all 
pixel signals) are calculated for each cluster. For the studies presented in this paper, we required the 
cluster energy to be >1 keV, which guarantees a negligible probability of accidental clusters from 
readout noise. Selection criteria specific to the different analyses will be described in Sections 0 
and0. 

Two sets of data were collected, with different operating parameters for the CCD readout. The 
DAMIC readout employs a correlated double-sampling technique [jl^ , where high-frequency noise 
is suppressed by measuring the voltage of the output node over relatively long intervals of time. The 
maximum pixel charge that can be recorded is limited by the dynamic range of the digitizer, since 
the pixel pedestal signal is proportional to the integration time. In the standard DAMIC readout, 
which is optimized for low noise, an integration time of 40 /f s per pixel corresponds to a maximum 
energy of ~20 keV. Radiogenic a particles have a range comparable to the pixel size, which leads 
to large localized charge deposits that can saturate the digitizer. Thus, a “low-gain” readout mode 
with an integration time of 0.6 /f s per pixel was used for runs dedicated to a spectroscopy. Due to 
the shorter integration time, this set of data presents a higher pixel noise (~40eV RMS), which is 
still negligible for MeV a particles. 

Table |T] summarizes the data collected for this analysis. We used the low-gain data (0.6 /is 
integration time) to evaluate the and ^^^Th contamination in the CCD by a spectroscopy 
(Section ||). The rest of the data (40/is integration time) was taken after the installation of the 
ancient lead shield, and presents a tenfold decrease in the background of electron-like tracks. The 
rate of a particles did not change, since the CCDs and surrounding components inside the box, 
where the as must originate, were not touched. We used this data set to constrain the activities of 
^^Si and in the bulk through a search of spatially correlated pairs of /3 tracks (Section 

3. Limits on uranium and thorium contamination from a spectroscopy 

Uranium and thorium and their decay products (Figure H) are ubiquitous in nature. They may be 
found within the bulk silicon of the CCD or the electrical elements of the device. They may also 
be present on the surface of the CCD, in particulates (e.g. dust) or as adsorbed radon daughters 
(mainly via deposition of ^*^Po from ^^^Rn decays in the air surrounding the device). 

The CCD is >99.5% electronic-grade silicon [1T311 by mass. The remainder consists of light 
elements used in the oxidation, metallization and doping of the silicon to fabricate the semicon¬ 
ductor structures necessary for operation. Due to the extremely high chemical purity of the silicon. 
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Figure 4. and decay chains. Alpha (j3) emitters are colored yellow (teal). For each isotope, a 
energies, j3 Q-values and half-lives are given. Isotope sequences expected to be in secular equilibrium are 
grouped by dashed lines. 


the substrate (photosensitive volume in Figure |T]) can be doped with a donor density as low as 
lO'^ [^. Thus, even an unexpectedly large contamination of or ^^^Th with an atomic 
abundance ~10% of the donor density would only correspond to ppt (10^'^) levels by mass. 

The CCDs are consistently handled in ESD-safe clean rooms (class 1000 or better). DAMIC 
CCDs are packaged at Fermilab, with facilities developed for DECam, an astronomy camera whose 
requirements for the deposition of dust on the CCD surface are much more stringent than that of 
DAMIC. However, the clean rooms are not radon-free, and deposition of radon daughters, mostly 
on the front CCD surface, will occur during the packaging. In particular, after being glued to their 
support, CCDs are left for two days under an air column of a few centimeters in height to cure 
the epoxy. We estimate that this procedure may induce a residual surface activity of and 
its daughters of less than 3x10^^ cm^^d^^ assuming a typical ^^^Rn activity for indoor air of 
30Bqm 

Many of the isotopes in the and ^^^Th decay chains (Eigure ^ are a emitters, and can be 
efficiently identified by a spectroscopy. 
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(a) Plasma a 


(b) Bloomed a 


Figure 5. a) A highly-diffase, round cluster due to plasma effect from an a particle originating in the bulk 
or the back of the CCD. b) An example of blooming, where an a particle originating in the front of the CCD 
produces a mostly vertical cluster. 


3.1 Characteristics and selection of a-induced clusters 

Radiogenic as lose most of their energy by ionization, creating a dense column of electron-hole 
pairs that satisfy the plasma condition [ITTll . The local electric field within the plasma is much 
greater than the electric field applied across fhe subsfrafe. For inferacfions deep in fhe subsfrafe, 
where fhe elecfric field is only along z, fhe charge carriers diffuse laferally foward regions of lower 
charge concenfrafion unfil fhe subsfrafe elecfric field becomes dominanf. Thus, as originating 
in fhe bulk or fhe back surface of fhe CCD lead fo highly-diffuse, round clusfers of hundreds of 
micromefers in diamefer. On fhe ofher hand, a parficles fhaf sfrike fhe fronf of fhe CCD deposif 
fheir energy less fhan 20 /im below fhe gafes. The high-densify charge cloud can easily overcome 
— and spillover — fhe barrier phases befween verfical pixels, while if is harder for if fo overcome 
fhe pofenfial barrier of fhe verfical channel slops befween columns. This phenomenon is known 
as blooming []T21], and leads fo mosfly vertical clusters. Examples of “plasma” and “bloomed” as 
detected in DAMIC are shown in Figs. |5(a)| - |5(b)| . 

Due fo fheir distinctive fealures, simple criteria are sufficienl fo efficienlly selecf and classify 
as. In fhe energy range of radiogenic as (>1 MeV), eleclron fracks are long (~mm or more) and 
deposif fheir energy in extended “worm”-like fracks over many CCD pixels. To differenliale elec- 
Irons from as, we determine fhe smallesl reclangular box lhal can conlain a cluster, and compute 
fhe fraction of pixels, fp^, in fhis “bound box” which are parf of fhe cluster. For small, symmef- 
ric clusters (i.e. as) fp^ is large (~7r/4 for a round cluster). For fhe long and irregularly shaped 
worms characferisfic of elecfrons, fp^ is small and decreases wifh increasing eleclron energy. Fig¬ 
ure |6(a)| shows fhe successful separalion befween j3s and as in fhe low-gain dala sel according fo 
fhis variable. 

To furfher dislinguish plasma from bloomed as, we calculate fhe spatial RMS Ox^y of fhe 
clusfer as fhe signal-weighfed RMS value of fhe pixels’ x,y coordinafe. Plasma as presenf a round¬ 
shaped cluster, wifh OxjOy ~ 1, while bloomed as are generally longer along fhe y axis, giving 
Ox/Oy < \.\n addition, fhe diffused clusfers from plasma as have more pixels (Apix) fhan bloomed 
as. We used fhe variable Apix Ox/Oy fo separate plasma from bloomed as (Figure |6(b)[). 


















(b) Plasma-bloomed a selection 


Figure 6. Selection of a particles, a) The fraction /pix as a function of the cluster energy. Clusters in the 
region above the dashed line are selected as as. b) The variable A^pixCTx/o’v as a function of the a energy. 
Plasma (bloomed) as are indicated by red (black) dots above (below) the dashed line. Open (closed) dots 
refer to clusters detected in the 675 fim (500 /rm) -thick CCD. Clusters are more diffuse in the thicker CCD, 
resulting in larger A^pix. 


Table 2. CCD physical properties and rate of observed as for the three CCDs installed at SNOLAB. 


CCD 

Mass / g 

Area / cm^ 

Bloomed rate / d ^ 

Plasma rate / d ^ 

500 /xm top 

2.2 

19 

0.87±0.17 

0.21±0.09 

500 /I m bottom 

2.2 

19 

0.87±0.17 

0.14±0.07 

675 /im 

2.9 

19 

0.63±0.15 

0.14±0.07 

Average 

2.4 

19 

0.79±0.10 

0.16±0.04 


3.2 Limits on radioactive contamination from a analysis 

The measured rate of as in the three installed CCDs is reported in Table 0. 

Most of the bloomed as are clustered around the characteristic energy of ^'°Po decay (5.3 MeV, 
Figure |6(b)|) , which may be present as residual surface contamination following exposure to ^^^Rn. 
Also, a significant number of as have energies <4 MeV, lower than any as from the and 
^^^Th chains (Figure Most likely, these are a particles which lose some energy before reaching 
the active region of the device, and originate from surface contamination of the CCD or nearby 
materials. If we conservatively assume that all bloomed as with energies <6 MeV are due to 
^'^^Po decays from ^'^^Pb contamination on the front surface of the CCD, we obtain an activity of 
0.078±0.010cm^^ d^'. This is twenty times larger than the upper limit of 3x 10^^ cm^^ d^^ pre¬ 
viously estimated from the exposure to ^^^Rn during packaging. However, a significant number 
of the observed as are likely to originate from the surfaces facing the CCD, i.e. the copper plates 
above the 675 /im and the top 500 /im CCD, and the silicon support piece above the bottom 500 /im 
CCD. Unfortunately, contributions from all these different surfaces cannot be disentangled with the 
available data. Similar considerations can be applied to plasma as, which should release their full 
energy if occurring in the bulk. Many of the plasma as have energies lower than those from 
and ^^^Th chains, suggesting surface contamination. Taking again the conservative assumption that 
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all plasma as with energies <6MeV originate from ^^°Po back-surface contamination, we obtain 
an activity of 0.012±0.004cm^^ 

Spectroscopy of plasma as can be used to establish limits on ^'°Pb , and ^^^Th contam¬ 
ination in the bulk of the CCD. Four plasma as whose energies are consistent with ^^‘^Po were 
observed (Figure |6(b) ). One of them cannot be ^'^’Po, as it coincides spatially with two higher 
energy as recorded in different CCD exposures, and is therefore likely part of a decay sequence 
(Section ^3]). When interpreting the other three as bulk contamination of ^'°Po (or ^'‘^Pb), an up¬ 
per limit of <37 kg^' d^' (95% CL) is derived. In the chain, the isotopes ^^°Th and 
^^^Ra decay by emission of as with energies 4.7^.8 MeV (Figure Since the isotopes’ lifetimes 
are much longer than the CCD exposure time, their decays are expected to be uncorrelated. No 
plasma as were observed in the 4.5-5.0MeV energy range, and an upper limit on the con¬ 
tamination of <5 kg^^ d^* or <4ppt (95% CL) is correspondingly derived (secular equilibrium of 
the isotopes with was assumed). In the ^^^Th chain, the timescale of the short-lived decay 
sequence of ^^^Ra-^^'^Ra-^^^Po is ~1 min, which is much smaller than the CCD exposure time. 
Thus, such a decay sequence would result in a single cluster when occurTing in the bulk, with a 
total of 18.8 MeV of energy deposited by the pile-up of the thr'ee decays. No cluster with energy 
>16 MeV is observed, which results in an upper limit of <15 kg^* d^^ or <43ppt (95% CL) on 
^^^Th contamination in the CCD bulk. 


3.3 Observation of spatially correlated a decay sequences 

We detected four plasma as with energies >5.5 MeV (Figure |6(b)| ), which cannot be accounted for 
by ^^‘^Po decay. They all occur in the top 500-/im CCD. Two of them, along with an a close in 
energy to the ^^‘^Po line, have their centroids within 1 pixel on the CCD x-y plane (Figure |7(a)| ). 
Since the accidental probability of this occurence is negligible, these three strongly spatially cor¬ 
related events must have a common origin. A likely explanation is that we have observed a decay 
sequence starting with a ^^^Th nucleus located in the thin (60 nm) ITO layer covering the backside 
of the CCD. The energies of the as and the time separation between decays are consistent with 
those from the decays of ^^^Th, ^'^Po and ^*^Po (Figure The two other as of the sequence, 
^^■^Ra and ^^°Rn, must have been emitted away from the CCD and thus gone undetected. Note 
that the recoiling nuclei corresponding to the ^^"^Ra and ^^^Rn decays may penetrate the device, 
but their deposited energy would be too small to be observable on top of the ^^^Po a once the 
detector resolution is considered. Likewise, the j3s from ^'^Pb and ^^^Bi must have been emitted 
away from the CCD or have deposited < 100-200 keV in the CCD. The observation of this single 
decay sequence is compatible with ~100ppb of ^^^Th contamination in the ITO [|T^. 

The highest energy plasma a (15.3 MeV) must be due to pile-up, either of two as (e.g. ^'^Po 
and ^^^Po) or as and j3s with energies < 100-200 keV (so that the beta track is too short to protrude 
outside of the a cluster). Likewise, the remaining 6.2 MeV a could be due to either a single a 
(e.g. ^^‘^Ra) or a less energetic a piled-up with a < 100-200 keV j3. In either case, these decays are 
unlikely to have occurred in the bulk of the CCD, as there are no isolated decays or short sequences 
of decays (<hours) corresponding to these energies in the or ^^^Th decay chains (Figure 

A spatial coincidence search between as and jSs is limited by the large number of accidentals 
from background, as the j8 rate is ~100 times greater than the a rate in the low-gain data set. 
In one of the latter data sets we have observed a coincidence between an a and a j3 within the 
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E = 5.39 MeV 


E = 6.75 MeV 


E = 8.66 MeV 




At = 17.8 d At<5.5h 

(a) Triple a sequence 



(b) a-jS coincidence 


Figure 7. a) Three a particles detected in different CCD images at the same x-y position. Their energies and 
the time separation between images are consistent with a sequence from a ^^^Th decay chain, b) A peculiar 
cluster found in a single image, consistent with a plasma a and a j3 track originating from the same CCD 
position. This may happen for a radioactive decay sequence occurring within the 8-hour exposure time of 
an image. 


same exposure, shown in Figure [7(b)| . As the energy of the a cannot be measured due to digitizer 
saturation, this particular decay sequence cannot be identified. 


4. Limits on ^^Si and contamination from /3 decay sequences 


210 


We have performed a search for decay sequences of two jS tracks to identify radioactive contam¬ 
ination from ^^Si and ^^°Pb and their daughters, whose p spectra extend to the lowest energies 
and could represent a significant background in the region of interest for the WIMP search. These 
isotopes do not emit a or penetrating /radiation, and their decay rates are significant for extremely 
low atomic abundances due to their 10-100 y half-lives, making conventional screening methods 
ineffective in determining their presence at the low levels necessary for a WIMP search. 

^^Si is produced by cosmic ray spallation of argon in the atmosphere, and then transported to 
the Earth’s surface, mainly by rain and snow. Detector-grade silicon is obtained through a chemical 
process starting from natural silica. Therefore, the ^^Si content of a silicon detector should be close 
to its natural abundance in the raw silica. Spectral measurements of radioactive background in 


silicon detectors suggest a rate of ^^Si at the level of hundreds of decays per kg day [1^]. ^^Si leads 
to the following decay sequence: 


^^Si — >^^P + j3 with Ti/2 = 150y, Q — value = 227keV (4.1) 

32p—Q-value = 1.71 MeV (4.2) 

^'‘^Pb is a member of the decay chain (Figure ^ and is often found out of secular equilib¬ 
rium, as chemical processes in the manufacture of materials separate it from other daughters. 
It may also remain as a long-term surface contaminant following exposure to environmental ^^^Rn 
(Section H). ^^°Pb leads to the following decay sequence: 
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Figure 8. Candidate j3 decay sequence found in data. The first cluster was detected in an image taken on 
2014/08/05 and deposited 114.5keV of energy. A second cluster, with energy 328.0 keV, was observed in 
an image taken 35 days later. Both tracks appear to originate from the same point (yellow star) in the CCD 
x-y plane. 


2'Opb ^ ^ 22y, Q-value = 63.5keV (4.3) 

—)> + j3^ with Ti /2 = 5d, Q —value = 1.16 MeV (4.4) 

The ^'°Ph nucleus decays 84% of the time into an excited state of ^'°Bi which promptly 
releases its 46.5 keV of energy hy internal conversion of an atomic electron (1C) in 80% of the 
decays or hy emission of a y-ray in 4% of the decays. ^'°Po is itself radioactive and decays hy a 
emission. The possible contamination from ^^‘^Po in the CCD has been discussed in Section 

The intermediate nuclei, ^^P and ^’°Bi, are expected to remain in the same lattice site as their 
parent nuclei and throughout their lifetimes. Therefore, the j3s produced by each decay pair should 
originate from the same pixel (out of 8x10®) on the x-y plane of the CCD. Through a search for 
electron-like tracks starting from the same spatial position, individual ^^Si-^^P and ^'®Pb-^'°Bi 
decay sequences can be selected with high efficiency. We performed this search with the lowest 
background data set (Table |I]) in the 675 /fm CCD. Given the background level (~ 10 electrons per 
day in a CCD), the number of accidental coincidences among uncorrelated tracks are small for 
periods of time comparable to the half-lives of ^^P and ^^°Bi. A candidate decay sequence found 
in the data is shown in Figure to illustrate the search strategy. 

4.1 Search procedures for spatially correlated j8 decay sequences 

The first step in the search for decay sequences is to find fhe end-poinfs of fhe j3 hacks. The 
procedure is illusfrafed in Figure ^ Firsf, we find fhe pixel wifh fhe maximum signal in fhe clusfer, 
and we use if as a seed poinf. Then, for every pixel of fhe clusfer we compufe fhe lengfh of fhe 
shorfesf pafh fo the seed point, where the path is taken only along pixels that are included in the 
cluster. We refer to this as the “distance” from the seed point. The pixel with the greatest distance 
is taken as the first end-point of the track. Finally, we recompute the distance of every pixel from 
the first end-point, and take the pixel with the largest distance as the second end-point of the cluster. 

To find a jS decay sequence, we calculafe fhe distance from the end-points of every jS cluster 
in an image to the end-points of every j3 cluster in later images. Thus, for every pair of clusters we 
have four distances corresponding to each end-point combination. The minimum of these distances 
is defined as the “cluster distance.” The pair is considered a candidate for a decay sequence if the 
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Figure 9. Algorithm to find the end-points of a cluster, a) The pixel with maximum signal is chosen as seed 
point, b) The distance of each pixel in the cluster to the seed point (star) is computed, c) The pixel with the 
largest distance is chosen as the first end-point (star). Distances to the first end-point are calculated, and the 
pixel with the largest distance is taken as the second end-point (the reddest pixel). 


cluster distance is smaller than 20 pixels and the clusters have at least one pixel in common. We 
refer to the cluster in the earlier (later) image as the “first” (“second”) cluster. 

To reduce the number of accidental pairs, we impose additional criteria on the energy of the 
clusters and their time separation. For the sequence search, we require the energy of the 

first cluster to be <230 keV and the energy of the second cluster to be <1.8 MeV. For the 
^'*^61 sequence search, we require the energy of the first cluster to be in the range 30-65 keV, 
which mostly includes the 80% of decays with an electron from internal conversion (Eq. EH). 
Clusters from IC will be diffusion limited, as they will be constituted by a cascade of j3, conversion 
and Auger electrons with energies <30keV. The energy of the second cluster is restricted to be 
<1.2 MeV. 

Lastly, we require the time separation between the clusters of each pair, Atpair, to be less than 
five half-lives of fhe daughter nuclei. This corresponds fo 70 (25) days for fhe ^^Si-^^P (^^°Pb- 
^'^^Bi) decay sequence search. 

4.2 Pair selection efficiency 

The efficiency of the pair selection described in Section EH was estimated by Monte Carlo sim¬ 
ulations. We used the MCNPX5 [ jT^ program with full electron tracking to simulate j3 particle 
interactions in a rectangular block of silicon with the same dimensions of a DAMIC CCD. The 
energy deposited by a particle is recorded in a mesh with cells of the size of a CCD pixel. The 
average z-coordinate of the simulated particle track within each cell is also recorded. Then, a re¬ 
alistic resolution and diffusion model is applied to the number of ionized charge carriers from the 
position at which they were produced along the electron’s path, and a simulated pixel cluster on the 
x-y plane of the CCD is obtained. For the purpose of the decay sequence search, the origin of the 
j3 particle was randomly generated in the silicon volume of the CCD to simulate a decay occurring 
in the bulk. 

To properly include the measured readout noise in the analysis, we used data “blanks.” These 
are zero-length exposures read out immediately after every data exposure, which feature true read¬ 
out noise patterns but no physical tracks. There are 159 blanks taken during the decay sequence 
search period. In each of them, we introduced the simulated clusters of five jSs from ^^Si fo ap- 
proximafely reproduce fhe rate of electron-like tracks measured in data. For each simulated ^^Si 


























cluster distance / pix 


(a) Simulated sample search 



(b) eti„e for 32si-32p 


Figure 10. a) Cluster distance distribution (black histogram) for the candidate pairs found in the 

simulated sample. Accidental pairs are represented by the red shaded histogram, b) Probability for the ^^P 
decay to be identified within the real-time of the data taking period, as a function of the time of decay of the 
parent ^^Si. The probability is zero when data was not acquired. 

decay, a decay was generated from the same location in a later image, distributed in time ac¬ 
cording to the half-life of With this method, 496 decay pairs were introduced in the sample of 
blank images. We used the standard CCD image reduction (Section P3| ) to reconstruct clusters. To 
properly account for inefficiencies of the CCD, bad pixels found in data were also masked in the 
blanks. Then, the pair selection procedure was applied to this set of simulated images. We found 
504 candidate pairs (Apairs) with a cluster distance distribution shown in Figure |10(a)| . As most of 
the pairs have a cluster distance <7 pixels, we adopt this additional criterion for the pair selection. 

First, we estimated the efficiency Cpair of requiring a cluster distance <7 pixels and cluster 
energies consistent with a ^^Si -^^P decay sequence (Section P~T|) . Of the 486 pairs with cluster 
distance <7 pixels, 475 are true pairs. Thus, the pair selection procedure is highly efficient in 
recovering decay sequences (£pair=475/496=95.8%), while keeping background at a reasonable 
level. The inefficiency is mostly due to pairs for which one or both tracks were not properly 
reconstructed due to masked pixels (16 out of 21 lost pairs). 

Second, we determined the efficiency Ctime to selecf a ^^Si-^^P decay sequence with pairs 
separated in time by less than 70 days (Section P~T1 ) when the parent ^^Si decays during the live¬ 
time of our data set. In Figure |10(b)| , the probability for a ^^P decay to be identified within the 
live-time of the data taking period is calculated as a function of the time of decay of the parent ^^Si. 
Integrating over the live-time of the data set gives an average £time =51.3%. 

An analogous simulation study was performed for the ^^®Pb-^^°Bi decay sequence. In this 
case, £pair is mainly determined by the requirement for the energy of the first cluster to be in 
the range 30-65 keV (Section ^~T|) . To estimate this efficiency, we have considered all possible 
processes following ^*°Pb decay, including the cases where some of the energy is radiated in 7 - 
rays or X-rays that escape the decay site and will not form part of the cluster. Only 82% of ^'^^Pb 
decays fall in this energy range. Of those, about 7.4% are not selected because their measured 
cluster energy is <30keV due to pixel saturation. The corresponding pair selection efficiency is 
found to be £pair = 71.0%. Also, an average £tinie =65.1% is obtained for this decay sequence. 
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Figure 11. a) Cluster distance distribution for the 16 candidate pairs found in 56.8 days of data, 

b) Distribution of the time separation of the 13 pairs with cluster distance < 7 pixels. The red line is 
the expectation for accidental pairs. 




Figure 12. Distribution of the number of accidental pairs with cluster distance <7 pixels for the 
(a) and ^'‘^Pb (b) decay sequences. The distributions were obtained by performing 1000 searches over 

spatially randomized data images. The red-dashed line is the Poisson distribution corresponding to the mean 
of the histogram. The red entry is the number of candidate pairs found in data. 

The overall efficiency for detection of ^^Si -^^°Bi) decay sequences in our data set, 

£x=£pair£time> where X = Si and Pb, is determined to be esi=49.2% (£pb=46.2%). 

4.3 Limits on ^^Si and contamination 

The search procedure described in Section o was applied to the data. No candidate was found 
in the ^'®Pb-^^°Bi search. The ^^Si-^^P decay sequence criteria yielded 16 candidate pairs, 13 
of which have a cluster distance <7 pixels (Figure |1 l(a)[ ). The Afpair distribution of these 13 pairs 
is shown in Figure |1 l(b)| . Also shown are the corresponding distributions expected for purely 
accidental pairs, which were estimated directly from the data. For this purpose, we randomized the 
position of the clusters in all data images, effectively eliminating any spatial correlation between jS s 
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Table 3. Summary of results presented in this paper. All values are 95% CL upper limits or intervals, except 
for the ^*‘'Po surface rate, where the uncertainties are 1-ff. The two measurements of the ^*‘'Po surface rate 
correspond to the two (back, front) CCD surfaces. For and ^^^Th, we quote the corresponding ppt 
contamination in parentheses. 


Analysis 

Isofope(s) 

Tracer 

Bulk rafe 

Surface rale 

mefhod 


for 

kg^^ d^^ 

cm^2 j-i 

a 

210po 

210pb 

<37 

0.011±0.004, 0.078±0.010 

specfroscopy 

234u 230226 r^ 

238u 

<5 (4ppf) 

- 


224Ra-220Ra-216po 

232Th 

<15 (43ppl) 

- 

j3 spafial 

32si_32p 

32Si 

80;“^^ 

- 

coincidence 

210pb_210gi 

210pb 

<33 

- 


from decay sequences potentially present in the data. In the randomization process, the overall x-y 
distribution of the clusters was maintained to avoid introducing a bias. Accidental pairs were then 
found in the set of randomized images by applying the search criteria. The procedure was repeated 
one thousand times, obtaining for each randomized search the accidental pairs’ distribution of 
cluster distance and Atpair- The average of the one thousand distributions (red line) is shown in 
Figs. |ll(a)| - |nX51 . Within the limited statistics, an excess of pairs at small cluster distances is 
observed in data. Their Atpaii distribution is compatible with the 14-day half-life of 

To evaluate the significance of the excess in the search and establish a limit on ^'^^Pb 

contamination, the probability distribution of the number of accidental pairs, Pacc^ was determined 
with the same spatial randomization procedure described above. For each randomized search, the 
number of accidental pairs with cluster distance <7 pixels, Aacc, was obtained. A Aacc distribution 
was derived from one thousand randomized searches, and taken as an estimate of Pace- We verified 
wifh Monfe Carlo simulafions fhaf fhe hue Pace is indeed recovered by fhis procedure. Figure |T^ 
shows fhe Aacc disfribufions obfained for fhe ^^Si-^^P and ^^°Pb-^^°Bi decay sequence searches. 
These disfribufions are well represenfed by Poisson disfribufions wifh mean of 6.5 for fhe ^^Si- 
^^P search, and mean of 0.8 for fhe ^'^^Pb-^'^Bi search. Wifh fhis assumption, fhe 13 observed 
candidate pairs correspond fo 1.2<Asi<15.3 (95% Cl) uni], where Asi is fhe number of esfimafed 
^^Si-^^P decay sequences in fhe dafa. The null resulf for fhe ^^°Pb-^^°Bi search corresponds fo 
Apb<2.5 (95% CL) QTTl], where Apb is fhe number of estimated decay sequences in fhe dafa. 

The ^^Si and ^'‘^Pb decay rates are fhen obfained as Ax/£x/T/Mccd! where ex is given in 
Section T=56.8 d is fhe dafa live-time and Mccd=2.9 x 10^^ kg. We estimate a decay rate of 
80lg5°kg^' d^^ (95% Cl) for ^^Si in fhe CCD bulk. This resulf also esfablishes fhe defection of 
^^Si -^^P pairs in our dafa sef af 98% CL. The derived upper limif on fhe ^'®Pb decay rafe in fhe 
CCD bulk is <33kg-‘ d-i (95% CL). 

5. Conclusions and Outlook 

We have presenfed novel analysis mefhods fo measure radioaefive confaminafion in fhe high- 
resisfivify silicon CCDs of fhe DAMIC experimenf. We exploifed fhe unique signafures of a 
particles in CCDs fo perform a specfroscopy and search for nuclides of fhe uranium and fhorium 
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chains. Also, we searched for and in the bulk silicon of the CCD by identifying pairs of 
spatially correlated j8 tracks compatible with and ^'°Pb-^^°Bi decay sequences. 

The results are summarized in Table We placed stringent 95% CL upper limits on the 
presence of radioactive contaminants in the silicon bulk. The and ^^^Th decay rates were 
found to be <5 kg^‘ d^^ and <15 kg^* d^\ respectively. Also, we established an upper limit of 
~35 kg^' d^* (95% CL) on the ^'‘^Pb decay rate, obtained independently by a spectroscopy and the 
P decay sequence search. In addition, we have measured the decay rate of ^^Si in the silicon bulk 
to be kg^'d^^ (95% Cl). Since we detect single nuclear decays with high efficiency, our 

analysis methods have near-optimal sensitivity to radioactive contamination for a given exposure. 
In particular, the capability to identify ^^Si-^^P and decay sequences, unique among 

particle detectors, allows us to measure these contaminations at levels that are orders of magnitude 
lower than would be possible with any other screening technique. 

The rate of as on the surface of the CCDs was also measured, and results are reported in 
Table 0. Since the observed rates may be partly caused by contamination on surfaces surrounding 
the CCDs, these results should be taken as upper limits on the surface contamination of the CCDs. 

These levels of radioactive contamination will allow DAMIClOO to probe WIMP-nucleon 
spin-independent interaction cross-sections as small as 10^^ pb for WIMPs with masses as low as 
2GeV/c^. In addition, material screening and handling procedures implemented for DAMIClOO 
should result in a tenfold reduction of the present background level, allowing for even more strin¬ 
gent limits on the radioactive contamination. Confirming the measurement, currently statistically 
limited, of ^^Si in high-resistivity silicon will be particularly relevant. In fact, the presence of the 
low energy P decay of this cosmogenic isotope may impose additional constraints on the next gen¬ 
eration WIMP searches with high-purity silicon detectors, including the identification of a source of 
silicon with low ^^Si content to fabricate the detector and the requirement of ^^Si-^^P single-event 
identification for background suppression. 
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